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We study the optimized perturbation theory (OPT) at finite temperature, which is a self- 
consistent resummation method. Firstly, we generalize the idea of the OPT to optimize the 
coupling constant in \(f> 4 theory, and give a proof of the renormalizability of this generalized 
OPT. Secondly, the principle of minimal sensitivity and the criterion of the fastest apparent 
convergence, which are conditions to determine the optimal parameter values, are examined 
in X(f> 4 theory. Both conditions exhibit a second-order transition at finite temperature with 
critical exponent /3 = 0.5 in the two-loop approximation. 

(N ' 

. §1. Introduction 

Experiments on ultra-relativistic heavy ion collisions at RHIC (the Relativistic 
Heavy Ion Collider) and LHC (the Large Hadron Collider) are being planned f< 
the purpose of producing hot matter similar to that realized in the early universe. lJ 
In such a high-temperature (T) phase, quarks and gluons are deconfined and chiral 
symmetry is expected to be restored. The dynamics of such a phase transition from 
<j_) | hadronic matter to a quark-gluon plasma is described by quantum chromodynamics 

(QCD) at finite T. 

At present, there are two major calculational methods in QCD at finite T. One 
is the numerical simulations in lattice QCD.0) This method has revealed that the 
restoration of chiral symmetry takes place at T c ~ 150 MeV. For two flavors massless, 
the phase transition has been shown to be of second order. Analysis of the case with 
two light flavors andone medium-heavy quark, which corresponds to the real world, 
is also under study cP. One of the disadvantages of lattice QCD is the difficulty of 
treating? real-time excitations at finite T in a straightforward way on the Euclidean 
lattice, fy 

Another method of treating^ QCD for T ^ is based on the hard thermal loops 
resummation scheme (HTLRS).tP It is generally known that naive perturbation the- 
ory breaks down at high temperature,!])'!!' even if the coupling constant is small. 
This is the reason that we must resum higher-order terms at finite T. HTLRS is a 
method of resuming the higher-order terms at high T up to O(gT) (g is the QCD 
coupling constant) when the external momenta are q ~ O(gT), a situation referred to 
as 'soft'. With this method, physical quantities such as the soft gluon damping rate 
and the dilepton production rate have been calculated in the quark-gluon plasma 
(QGP) phase. B> Since it is an effective resummation method which works only at 
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high T, it cannot be applied to systems at low T. Also, in theories with spontaneous 
symmetry breaking (SSB), the loop expansion is more relevant than weak-coupling 
expansions, such as HTLRS. 

To treat non-perturbative QCD effects at low energies, related to spontaneous 
symmetry breaking, effective models, such as the linear a model, have been devel- 
oped. In such approaches, however, the serious problem of tachyonic poles appear©, 
when one studies symmetry restoration at finite T. This problem shows up even 
below T c , and it causes the breakdown of the thermal perturbation theory. There- 
fore, a resummation method applicable to a wide range of T from low T to high T 
is required when effective models are studied. 

To this time, several methods have been proposed for self-consistent resumma- 
tion. IHP In these methods, there appear two major problems: EiP (i) appearance of 
T-dependent ultraviolet divergences, which makes the renormalizability non-trivial 
at finite T, and (ii) the breakdown of the Nambu-Goldstone (NG) theorem at finite 
T when SSB occurs. 

In Ref. tlP, we have generalized the optimized perturbation theory (OPT) at fi- 
nite T, which is a resummation method for solving the above mentioned problems. 
A similar idea has been proposed under the names of the delta-expansion, the varia- 
tional perturbation theory, and so on. These ideas are applied to field theory at finite 
T in Ref.Ej) andli^. In Ref. t2P, the O(N) </> 4 theory with the principle of minimal 
sensitivity (PMS) (see §2.1) is studied. Since the formulation in that work is not 
based on the loop expansion, its authors could not prove renormalizability Also, the 
NG-theorem is not satisfied in a straightforward way. In Ref.Ej), A(/> 4 theory with 
the fastest apparent convergence (FAC) condition (see §2.1) at high T is studied. 
There, the renormalization up to the two-loop order is carried out in the symmetric 
phase. However, a rigorous definition of the method at higher orders and a proof of 
its renormalizability are not given. 

Our method has the following advantages over the previous methods. 

1. Renormalization can be carried out automatically at each order of the OPT. 

2. The Nambu-Goldstone theorem is satisfied in any given order of the OPT. 

In this paper, we further generalize the OPT studied in Ref.0 to optimize 
the coupling constant in the A0 4 theory. Also, we investigate whether the phase 
transition is correctly described by the OPT in A</> 4 theory, which is known to have 
a second order transition. We examine the principle of minimal sensitivity and the 
criterion of the fastest apparent convergence, which are conditions to determine the 
optimal parameter values, in the two-loop approximation. 

The organization of this paper is as follows. In §|2|, we generalize the OPT 
introduced in Ref.E-if to optimize the coupling constant in A</> 4 theory. The proof 
of renormalizability in this generalized OPT is also given. In §|3|, we investigate 
the phase transition in A</> 4 theory in the OPT. The principle of minimal sensitivity 
(PMS) and the criterion of the fastest apparent convergence (FAC) @ at the two-loop 
order are examined. In §^|, we give summary and concluding remarks. 
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§2. Formulation of the OPT with optimization of the coupling constant 

In A0 4 theory with spontaneous symmetry breaking (SSB), the naive perturba- 
tion theory has the following difficulties. 

• Higher-order terms give larger contributions than lower order terms (the exis- 
tence of the hard thermal loop). 

• Tachyonic poles appear when one studies the restoration of symmetries. 

In the following, we introduce the optimized perturbation theory (OPT) with both 
mass and coupling constant optimization. This theory allows us to avoid the above 
problems. 

2.1. Formulation of the OPT 

We reformulate the OPT to incorporate the optimization of the coupling con- 
stant. This extension can be performed in the same way in Ref.EJ. Let us start 
with the renormalized Lagrangian £ of A</> 4 theory: 

£(0;/, 2 ,A) = ip0) 2 -^V]-^ 4 

+lA(X)(d<t>) 2 - \b(\)^4> 2 - \C{\)^ + D(\)^ (2-1) 

Here we have explicitlv written the arguments fj? and A in £ for later use. Since we 
adopt the M S scheme t3 (which is a mass-independent renormalization scheme) with 
dimensional regularization, A, B, C and D are functions of A only. For simplicity, we 
omit the dimensional factor K 4_n that multiplies A. Here, k is the renormalization 
point and n is the number of dimensions. 

Step 1 (Definition of the 5 expansion) 

The loop-wise 5 expansion for the effective action is defined as 



r 1 f h / T r 

r[ip 2 } =5hln / [#] exp — / d A x C((f> + ip;[i 2 ,\) + J<f> 
J Oh Jo L 



(2-2) 



where J = —dr[<p]/d<p and Jq^ T d A x = dr J d 3 x. At T = 0, it is equivalent 
to the naive h expansion. EJ) However, it does not coincide with the h expansion for 
finite T because h is also contained in the upper limit of the integral. 

The counterterms are also expanded in powers of 5. Because the renormalization 
is performed at T = 0, the counterterms are the same as those in the naive loop 
expansion. 

Note that 5 is introduced only to define the order of the perturbation theory. At 
the end, we set 5 = 1. 

Step 2 (Splitting the mass and coupling) 

The mass and the coupling constant are decomposed as 

2 2 / 2 2\ 2 

^ = rnr - [m L - ^) = m z - x, ^.3) 
A = 9 -(9- A) = g-rj. 
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Namely, we add and subtract the mass term m 2 and the coupling constant g, and 
we define \ = m 2 — [i 2 and rj = g — A. Using the expression in Eq.(2-3) to substitute 
into Eq.(2-1), we can rewrite the Lagrangian as 

£0; m 2 ,x, g, v) = £(4>;m 2 - x,g - v) 

= I[W) 2 -mV]-|^ 4 + ^ 2 + > 4 

+ ±A(g - rj)(d^ 2 - X -B(g - 7?)(m 2 - 

-^C(g - t?)0 4 + D(g - ^(m 2 - X 2 ) 2 . (2-4) 

It is important that the identities in Eq.(2-3) are used not only in the standard 
mass and coupling terms but also in the counterterms □ to make the order by order 
renormalization possible in the OPT. 

To obtain a non-trivial loop expansion, we need to assign 5 as 

m 2 = 0(l), A = 0(1), X = 0(5), V = 0(5). (2-5) 

Thus, the tree-level mass becomes m 2 +gip 2 /2 instead of /i 2 + \(p 2 /2 in the symmetry 
broken phase. The order of 5 is increased by inserting the new vertex x4> 2 : /2 or ryc/> 4 /4!. 
The physical reason behind the relations in Eq.(2-5) is the fact that x an d f] reflect 
the effect of interactions. 

Since Eq.(2-4) is simply a rearrangement of the parameters, the effective action 
should not depend on the arbitrary parameters m 2 and g. However, since we cannot 
calculate all orders in actual calculations, the physical quantities depend on artificial 
parameters. Methods for determination of these parameters are given in the next 
step. 

Step 3 (Determination of m 2 and g) 

One can determine_the optimal parameter values of m 2 and g using methods 
proposed by Stevenson. 

(a) The principle of minimal sensitivity (PMS): 

^=0, ^^0. (2-6) 
om og 

Ol iP' l 1 ) represents a physical quantity calculated up to L-th (L'-th) order. 
Since m 2 and g are artificial parameters added by hand, the physical quantities 
should not depend on them. 

(b) The criterion of the fastest apparent convergence (FAC): 

O L -O L _ n = 0, O'v - 0' L ,- n , = 0, (2-7) 

where n (n') is chosen in the range l<n<L(l<n< L'). This condition 
requires that the perturbative corrections in O'l (P' l 1 ) should be as small as 
possible for a suitable value of m (g) . 

* This point was first made in Ref. 0. 
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These conditions are reduced to self-consistent gap equations. Therefore, OPT cor- 
responds to a generalization of the mean field approximation. 

An extension of the Lagrangian to the O(N) case is straightforward. Actually, 
one can easily see that the Nambu-Goldstone theorem is fulfilled in completely the 
same way as in Ref.lHP. 

In the next section, proof of the renormalizability in this extended OPT is given. 

2.2. Proof of the renormalizability in the OPT 

Consider a naive m-th loop order renormalized correction to the effective action, 
pir™) _ Then, the effective action up to n-th order can be written as 



m=0 



Here, each is written as the sum of r^ m \ which does not contain any countert- 
erms, and C^ m \ which includes more than one counterterms: 

r M =r (m) +c (m)_ (2 . 9) 

Also, rjj m) can be expanded in the renormalized mass [i and the renormalized 
coupling A around fi 2 = and A = 0: 

oo I 



^^EEfoW^+O- ( 2 - 10 ) 

i=0 j=0 



Here 



(m) = 1 d i & (m) 

ij ~ iij\d(v 2 yd\j 



^2=0^=0 



(2-11) 
(2-12) 



Since 4 is finite, + must also be finite. 

Next, we consider the OPT. In this method, /j, 2 and A are rewritten in Eq.(2-3). 
Then, the expressions in Eq.(2-3) are used to substituted into (2-10), yielding 



4 m) = E E(^ 2 " XT(9 - r,) j (r^ + C[f). (2-13) 

Since m 2 is 0(1), A is O(l), \ is 0(6) and rj is 0(6), the order of 6 increases if terms 
include \i which is 0(6) and/or rj, which is 0(6). Thus, by expanding the parameters 
(m 2 — xT an d (g — rj) 3 in Eq. (2- 13), terms of higher order in 6 are generated: 

r M = r M (jm) + r W( S (m+-l)) + . . . 

oo I 

= E El(™ 2 )V - {*x(™ 2 rV + Mm 2 ) 1 ^- 1 } + ■ ■ -](rjf> + cpm 

i=0 j=0 
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where 

r ^\s m+a+b ) = ,C a {-x) a {m 2 r a 3 C h (-r 1 ) b ^- b (rif + C^f). (2-15) 
Thus, up to the m-th 5 order correction, reads 

m 

4 m = E4 m " S) (n- (2-16) 

s=0 

Since the renormalization is carried out up m- th order, each (0 < I < m) is 
finite. Thus, -Tj> m is also finite. 

If a theory does not have a mass term, the renormalization of the OPT is defined 
as the limit /i 2 — > of Eq.(2-16). By using this definition, one can obtain the Debye 
mass even when the original Lagrangian does not have a mass term. 

§3. Phase transition in A(/> 4 theory 

3.1. Calculation of the effective potential 

In this section, we apply the OPT to A(/> 4 theory. Firstly, we calculate the 
ffective potential for finite T up to the 2-loop level using the real time formalism. 
The coefficients of the counterterms for Eq.(2T) at the 2-loop levelEP are 



A(X) 



J? 1_ 

(4vr) 4 24e' 



A 1 A 2 / 1 1 
A 2 3 A 3 f 9 3 



(4vr) 2 2e (4vr) 4 v 4e 2 2e 
1 1 A 1 
D(A) -"(4^4i" (i^si 2 "' (3 ' 1} 

where we adopt the MS scheme. (The factor n^~ n ^ multiplying Alias been dropped, 
as ahovej The^effective potential for (2-1) is calculated in Ref.ElP at T = and in 
Ref.E3),B) and© for T ^ 0. The result is (see Appendix 

V = V° + V 1 + V 2 , (3-2) 

V° = + (3-3) 

y^-^d-^H r^Hi-e-n, 0-4) 



(4tt) 2 4 v 2 k 2 ' Jo (2tt) 2 f3 
V*=±K? + *f S , + *fc„ (3-5) 

where M = fi 2 + \ip 2 /2, f3 = l/T and £7 = V*: 2 + M 2 . The definitions of K u S s and 
Cf are found in Appendix |A| (Eqs.(AT4), (AT2) and (A-15), respectively). 
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3.2. Application of the OPT to A</> 4 theory 

Let us apply the OPT to Eq.(2-1). This leads to 

£(cf>;m 2 ,x,9,v) = \m? - m 2 <t> 2 ] - |<^ 4 + ~ X «£ 2 + ^ 

+ ^A(g - V )(d0) 2 - X -B(g - ^(m 2 - 

-^C(g-r))<f> 4 + D(g-r})(m 2 - X 2 ) 2 , (3-6) 

where we have used Eq.(2-3). New diagrams for the one-point functions and the 
vacuum energy diagrams are given in Appendix |B| in Figs. 11 12, respectively. 
Up to 0(5 2 ), the effective potential Eq.(3-2) becomes 

V = V° + V S + V S \ (3-7) 

T/-0 2 2 i 9 4 to a\ 

V =-m ip + , (3-8) 



' /A ' |ifc 2 ln(l-e-^), (3-9) 



o (2vr) 2 /3 

V S2 = U + ^ + lK t) K t + ^f(S s + C f ). (3-10) 

In A0 4 theory at high T, only the tadpole diagram in Fig.l (A) is the hard 
thermal loop (HTL). Therefore, we must resum the cactus- type diagrams in Fig.l. 
Since they do not depend on external momentum, only the mass term is modified. 
Thus, the HTLs resummation is performed by shifting the mass term. This happens 
only when one considers a theory such as Xcfr model. If we consider QCD, all vertices 
with N external gluons and vertices with N — 2 external gluons and 2 external quarks 
identified HTLs which must be summed up. 



(A) (B) (C) 

Fig. 1. Tadpole and cactus diagrams in A0 4 -theory. In this theory, only a tadpole diagram become 
hard thermal loop. 

In the following subsections (3.3 and 3.4), we adopt the optimization only for 
the mass term in A0 4 theory for simplicity. We consider the case fi 2 < 0, and the 
restoration of the symmetry is discussed under various conditions (PMS and FAC, see 
§2.1). The full OPT case (optimization of the mass and coupling terms) is discussed 
in S3. 5. 
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3.3. The PMS condition 

Here, we investigate the principle of minimal sensitivity (PMS) in one-loop and 
two-loop orders. Since we study the static nature of the phase transition, the thermal 
effective potential V(tp, m 2 ) is chosen as the relevant physical quantity Ol in Step 3 
of §2.1. 

3.3.1. 1-loop analysis 

For the optimal condition in Step 3 of §2.1, we adopt the following condition at 
the 0{8) level: 

dV 0+5 (ip,m 2 ) 

where V 0+S ((p, m 2 ) = V°+V s (see Eq.(B-l) in Appendix [B]) . However, this condition 
does not lead to the appropriate gap equation. 

Differentiation with respect to m 2 corresponds to cutting one of the internal 
lines of V(<p, m 2 ). This is because the power of the propagator is raised by 1. As 
one can easily see from Fig.|6] in Appendix |A|, cutting the internal line of V((p,m 2 ) 
cannot produce HTLs (Fig.l (A)). Therefore, Eq.(3Tl) cannot sum the tadpole type 
diagrams, and it is not meaningful to adopt the PMS condition at the 0(5) level. 
Thus, we need to go to the next order, which is two loops. 

3.3.2. 2-loop analysis 

The PMS condition for the 2-loop effective potential reads 

om z 

where v 0+5+s2 (if, m?) = V° + V s + V s2 . Cutting one of the internal lines of Fig.|lO| (a) 



leads to HTLs in the scalar theory. The explicit form of (3-12) is given in Appendix 
B.l. At high T (in the symmetric phase), Eq.(3T2) reduces to 



dV°+ s+s2 (ip,m?) 



dm 2 

This equation gives the solution 



Sm 2 -^0 XT 2 T 

- k-s-W < 343 > 

<p=0 



AT 2 

m 2 (T)= M 2 + ^ r . (3-14) 

Since Eq.(3T4) corresponds to the Debye screening mass at high T, the condition 
(3-12) correctly resums higher order terms and recovers the reliability of the pertur- 
bation theory at finite T. 

To determine the vacuum, we must also solve the equation 

dV(y,m 2 ) 

^ = 0. (3-15) 

ip=ip 
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The derivative with respect to tp does not act on m 2 by definition, even if m 2 (ip) 
depends on ip. However, the gap equation Eq.(3T2) leads to 



Thus, in this case, the total derivative with respect to p is equal to the partial 
derivative. This is an advantage of the PMS condition over the FAC condition 
in studying the phase transition if m 2 (p), which is the solution of Eq.(3T2), has 
physical meaning for all p. 

Initial condition 

We solved Eqs.(3-12) and (3-15) numerically. There are three parameters in these 
equations, fi 2 , A and k. (m 2 is determined by Eq.(3-12).) Since we assume that the 
loop expansion at T = is a valid approximation, the renormalization point k is 
chosen so that fj? = m 2 is satisfied. Thus, the OPT is equivalent to the naive-loop 
expansion at T = 0. In other words, we use the OPT only for resummation at finite 
T to recover the reliability of the perturbation theory at finite T. (Note, however, 
that the relation fi 2 = m 2 at T = must be satisfied in Eq.(3T2).) Although the 
explicit values of the parameters are not important for our qualitative study (in fact, 
these parameters are normalized by their initial values), the initial values, A = 10.0 
and p = 10.0 were used for simplicity. As a result of solving Eqs.(3-12) and (3-15) 
simultaneously, we obtain /i 2 = m 2 = —170, k 2 = 87.6 and M 2 = m 2 + ^p^ = 330. 

Results of numerical calculation 

The results of the numerical calculations with the PMS condition are shown in 
Fig.2. Figure 2 (A) plots the tree-level mass M 2 (T) = m?(T) + \p$/2 with the 
left vertical scale and the optimized parameter m 2 (T) with the right vertical scale. 
M 2 (T) is clearly not tachyonic for any T. This result, Eq.(3-14), confirms that OPT 
with the two- loop PMS condition for V(ip,m 2 ) is successful for the resummation of 



Figure 2 (B) shows the T dependence of the thermal expectation value po(T) 
divided by ipo at T = 0. From this result, the phase transition is found to be of second 
order. Figure 3 (A) shows the second derivative of V(p,m 2 ) with respect to p at 
p = <po(T) as a function of T/T c . This figure also shows the second order nature 
of the transition. When the transition is of second order, the effective potential 
becomes flat for po(T) = at the critical temperature T c , which is expressed by the 
following equation: 




2 ) dm 2 
dp 



(3-16) 



HTLs. 



d 2 V(p,' 
dp 2 




ip=0 



= 0. 



(3-17) 



The effective potentials at T = and T c are shown in Fig.3 (B). We can also confirm 
the second-order nature of the phase transition from Fig.3 (B). We note here that 
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a second order_transition for the A</> 4 -theory is expected from the renormalization 
group analysis c3 and lattice simulations at finite T.E$ 

We found that the critical exponent (3, which is defined by 



<p (T) cx 



T-T r 



(3-18) 



becomes 0.5 in the two-loop analysis. This is the value expected from the Landau 
mean-field theory, which implies that our approximation is still within the level of 
the mean-field theory. Since the OPT corresponds to a generalized mean field theory, 
this was in fact anticipated. 

In Fig. 3 (C), the minimum value of the thermal effective potential V(ipo,m 2 ), 
which is equivalent to the Gibbs free energy, is shown. Its value decreases monoton- 
ically as T increases. 




Fig. 2. (A) The tree level mass M 2 (T) — m 2 + Xipo/2, with the left vertical scale and the mass 
parameter m 2 (T), with the right vertical scale obtained with the PMS condition. They are 
normalized by their vales at T = 0. (B) Vacuum expectation value <£o normalized by fo{T = 0). 
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3.4. The FAC condition 

In this subsection, we apply the FAC condition. The simplest condition to resum 
the HTLs in this case is 

E R {uj = Q,\k\ = 0;T) =0, (3-19) 
where Ur is a retarded two-point self-energy which is defined as 



ZR(u,k;T) 



(3-20) 



dip 2 

ip=tp 

In the following, we investigate the above condition at one- loop and two- loop orders. 

3.4.1. 1-loop analysis 

The HTLs can be resummed even with the one-loop FAC condition of Eq.(3T9). 
For this reason, this condition has been adopted in many studies for simplicity, and 
it allows the tachyon problem to be solved at low T. However, the phase transition 
becomes first order with this condition. Here, we do not recapitulate these results, 
because they are discussed extensively in the literature. Instead, we examine the 
two- loop condition Eq.(3T9) in §3.4.2 to compare the result with that for the PMS 
condition discussed in the §3.3.2. 

3.4.2. 2-loop analysis 

Since the physics should not depend on the artificial parameter m 2 , one expects 
that the result should not depend on the choice of the optimized conditions. However, 
the 2-loop PMS condition Eq.(3-12) leads to a second-order phase transition, and 
1-loop the FAC condition is known to lead to a first-order transition. Therefore, it is 
necessary to study whether the FAC condition at the two-loop level gives a second 
order phase transition. 

For the FAC condition at 0(S 2 ), we adopt 



S R (uj,k;T) = - 



dip 2 

¥>=¥>0 



0. (3-21) 



An explicit formula for V s2 is given in Appendix.^. At high T, Eq.(3-21) is reduced 
to 



d 2 v 52 



8<p 2 



pm^o XT . XT 2 . X 2 m 2 , , n 1, m 2 m 2 _ 

- i^ ( ^-ir ) + 2W {c+( " 2 + 2 ln ^ )ln ^ } 

95=0 v ' 

\ 2 T 2 T 2 

r (3.30-m-^). (3-22) 



24(4tt) 2 v k 2 

The first term on the right hand side (r.h.s.) produces the resummation of the HTLs. 
Initial conditions 

The initial parameters (at T = 0) are chosen as fi 2 = m 2 , A = 10.0 and ipo = 10.0 
as in the previous section. The condition fj? = m 2 is needed for agreement with the 
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naive loop-expansion at T = 0. Note that we recognize that the loop expansion at 
T = is a valid approximation. The other resultant values, which are determined 

-166, k 2 = 137 and 



by solving Eqs.(3T5) and (3-21) simultaneously, are fj, = m 
M 2 = m? + ±tpl = 334. 

Results of the numerical calculation 

The results are given in Figs. 4 and 5. One can see that the qualitative features 
of the figures are the same as the PMS results in Figs. 2 and 3 From Fig. 4 (A), 
one can see that the tachyon problem is also avoided in this case. Figures 4 (B), 5 
(A) and 5 (B) show the second order phase transition with (3 = 0.5. In this case, 
the Gibbs free energy decreases uniformly (Fig. 5 (C)). Thus, the 2-loop condition in 
both PMS and FAC give qualitatively the same results. This is the desired property, 
and it shows the validity of the OPT. 




Fig. 4. (A) The tree level mass M 2 (T) = m 2 + Xipo/2, with the left vertical scale, and the mass 
parameter m 2 (T), with the right vertical scale, in the case of the FAC condition. They are 
normalized by their values at T — 0. (B) Vacuum expectation value (po normalized by <po(T = 0). 




T I T c cp I q> {) { T=0) T I T t 

Fig. 5. (A) Second derivative of V(T) with respect to if for the FAC condition. (B) Effective 
potential at T = and T = T c . (C) Minimum value of the effective potential as a function of T. 



* The physical reason for the shoulder structure around T/T c ~ 0.7 in Figs. 4 (B) and 5 (A) is 
not yet understood. 
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3.5. Some remarks 

We tested the two conditions represented by Eqs.(3T2) and (3-21) in the OPT to 
study the phase transition for Xcj) theory. The two conditions give qualitatively the 
same results and show that the resummations were successfully performed. However, 
there are three remarks in order, which we give below. 

F AC with 0(6) + 0(5 2 ) 



One may try to use the two-loop FAC condition 

E s + s \u;,k;T) 



5+5*,.. d 2 V^ 2 (p,m 2 ) 



dp 2 



(3-23) 



instead of Eq.(3-21). This condition implies that all the loop corrections vanish 
at zero external momentum. In other words, m 2 contains all contributions up to 
two-loop order. From Eq.(3-23), the following equation is obtained: 

^_ = M W + ^. (3-24) 

At T = T c , if one assumes a second-order phase transition, Eq.(3-24) must be zero; 
that is m 2 (= M 2 ) = 0. However, the vanishing of the tree-level mass causes an in- 
frared divergence in the loop integrals. Actually, the left-hand side (l.h.s) of Eq.(3T5) 
diverges as M 2 — > 0. Thus, M 2 = is never satisfied at T = T c . From this argument, 
we see that a second-order phase transition cannot be realized with the condition 
given by Eq.(3-24). 

Full OPT 

The full OPT, which includes an optimization of the coupling constant, allows 
the possibility not only of avoiding the above infrared problem but also of going 
beyond the mean field approximation near T c . Therefore, it is very interesting to 
explore it. However, the PMS condition for V(p, m 2 ,g) requires tree- loop calculation. 
This can be understood as follows. Suppose one chooses the PMS condition at the 
two-loop level as 

dV + s + s \p,m 2 ,g) dV°+ s + s \<p,m\g) 

o 2 = > » = °' V 3 ' 25 ) 

In the symmetric phase {p> = 0), Eq.(3-25) can be reduced to 



dV 



dm 2 
dV_ 
~8g~ 



v=0 dM 2 



<>kt (x + 9Kt), (3-26) 



Kf. (3-27) 



first-order phase transition. (See also, Ref.Ej'.) 



* In Ref.e, Eq.(3-23) is used as the-EAC condition. Thus, the calculation there leads to a 
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From Eq.(3-27), Kt = is obtained. This gives a solution for m 2 . (Note that K 
depends on g only through M 2 = m 2 + gip 2 /2.) Then, this solution is substituted 
into Eq.(3-26), giving 



dV(<f,m 2 ,g) 



dm 2 



dM 2 

(p=0 



0K t 

X = 0. (3-1 



Since this relation also does not depend on g, we cannot determine g. Thus, Eq.(3-25) 
cannot give a solution for m 2 and g simultaneously in the symmetric phase. 
On the other hand, in the three-loop calculation, V(tp,m 2 , g) has the term 

This leads to the HTLs (tadpole diagram) in dV/dg, and depends on g even in 
the symmetric phase. Thus, the full OPT with the PMS condition for the tree-loop 
effective potential has the possibility of giving a solution for m 2 and g simultaneously. 

What about the FAC condition in full OPT? Suppose we use the following FAC 
condition: 

d 2 V s2 (cp,m 2 ,g) d*V s2 (<p,m 2 ,g) 

V " U ' [ ' 

Unfortunately, we could not find a solution which improves the previous results 
near T c at the two-loop level. We have tried all possible variations (V s2 (<p, m 2 ,g) 
in Eq.(3-30) replaced by V s+s2 (cp, m 2 , g), V s (tp,m 2 , g), and so on), but improved 
solutions were not obtained for the critical exponent. Thus, two-loop order, it seems 
that it is not possible to improve the previous analysis near T = T c in the full OPT 
with the FAC condition. In any case, further study along this line is necessary. 

Limiting temperature 

For sufficiently high T, there are no solutions with (3-12) and (3-21) for m 2 (T), 
because logarithmic terms of the form Iii(T/k) dominate. This means that the 
renormalization point k, which is fixed at T = 0, becomes a bad choice as T increases. 
To avoid this situation, one may try a renormalization group improvement. Since 
the typical "scale" of this system is T, one may choose k = T. In this case, no 
large m(T//c) appears. This renormalization group improvement would extend the 
applicability of OPT. However, in non-asymptotically free theories, such as X(p 4 
and 0(N) Xcp 4 theories, there eventually appears a Landau pole where the running 
coupling constant X(k = T) diverges at a certain T.EZP'ta' Thus, the theory is not 
applicable beyond this T. 

§4. Summary 

In this paper, we have generalized the optimized perturbation theory (OPT) 
at finite temperature (T) in the A</> 4 theory to incorporate the optimization of the 
coupling constant. The naive loop expansion for theories with spontaneous symmetry 
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breaking is known to break down at low T (tachyon pole problem) and high T 
(hard thermal loops (HTLs) problem). By contrast, we have shown that the OPT 
does not suffer from such difficulties. This is because that the OPT resums the 
higher-order terms of the naive perturbation theory in a consistent way by imposing 
appropriate conditions, such as the principle of minimal sensitivity (PMS) and the 
fastest apparent convergence (FAC). 

The advantages of the OPT over other self-consistent methods are that one can 
carry out the renormalization of ultraviolet divergences systematically and that the 
Nambu-Goldstone (NG) theorem at finite T is trivially satisfied at any given order. 

In §[|, we applied the OPT to Xcj) 4 theory to examine whether it can describe the 
finite T phase transition correctly. Carrying out the two-loop computation of the 
effective potential in the OPT, we have found that both PMS and FAC conditions 
give the correct second order transition. The critical exponent 0, however, is found 
to coincide with the mean-field value at this level. The full OPT, where both the 
mass and the coupling constant are shifted, may or may not improve the result. This 
remains as an open problem. 

Application of the idea of the OPT to gauge theories at finite T is also an inter- 
esting problem. However, to go beyond the hard thermal loops (HTLs) resummation 
scheme, one must solve the following two problems. 

• Since an infinite number of TV-point vertex functions (iV-gauge boson vertices 
and (N — 2)-gauge boson 2-fermion vertices) identified with HTLs in gauge 
theories and a naive mass term breaks the gauge symmetry, one must take 
into account an infinite number of effective vertices like HTLs resummation 
scheme. However, this infinite number vertices may create a difficulty for the 
renormalization. 

• It is known that HTLs resummation scheme breaks down at high T due to the 
so-called magnetic mass problem. EHP 

Therefore, we may need further generalization of the OPT to apply it for gauge 
theories. 
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Appendix A 

2-Loop Effective Potential in Xcp 4 -Theory 



g3re, we calculate the physical effective potential in (2T) in the real-time formal- 
> which is equivalent to the effective potential defined in the imaginary-time 
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formalism. Following Ref.li3, we compute the tadpole functions V^'°\ and then 
the effective potential is integrated over the classical field tp. Since the propagator 
in the T = part and part decouples in the real time formalism, calculations 

of the T = parts yield the same results as the ordinary T = perturbation theory. 
Hence, we can calculate the T-dependent term and T-independent term separately. 

The Feynman diagrams contributing to the pi tadpole functions up to the two- 
loop level are shown in Figs. 7, ^. The vertex with the number 1 (2) represents the 
type 1 (2) field self-interaction. The type 2 field is necessary to cancel pathological 
pinch singularities. 

The formulas of the tadpole functions read 

F (i,o) = y (l,0)l + F (i,o) 2; (Aa) 
- iV^ 1 = -i^V e f ^-iDf(k) + counterterms (A-2) 

2 J (2TT) n p 

V ilfi)2 = vtZ + V c ^ + V s ^, (A-3) 

r\ n k^ C rl n h n 



v (l,0)2 _ (~»A)V 4£ 
LV 2bubble ~ ■ K 



4 



+ counterterms, (A-4) 



2 



2 V ^ ^ ' J (2tt)' 

-^ 2 i^m (a-5) 

1,0)3 = t^J-J ^^[{iDiH^wf^iHh + k 2 ) 

-2iDf{ki)iDf{ki)iDf{k 2 )Wf{ki + k 2 ) 



+{iDf{k 1 )} 2 iDf{k 2 )iDf{ki + k 2 )\ 

(-i\(p)(-i\) 4e f d n h d n k 2 r u,, Vn ii,, Vn ii/ Z , , n 
+ g « y ( 2?r )n ( 27 r)n [t ^ ( l} ^ ( 2) * ^ ( 1 + 2) 



iDl 2 (ki)iDl 2 (k 2 )iDl 2 (ki + fe 2 )] + counterterms, (A-6) 



where, 



iDf (*) = ( ^ ; _ (A-7) 

1 _ e -/3|fcol/ 2 



+2vr5(A; 2 - u 2 



e /3 1 ^0 1 — 1 l_ e -^|fcol 
_ e -/3|fc |/2 1 
l_ e -/3|fc l e /3|fcol— 1 



* The "(1,0)" in V*- 1 ' ' indicates the 951 tadpole function. (pi is a type 1 field in the real-time 
formalism. 
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M 2 = ^ 2 + ^ 2 , n B {E) = (e? E -l)- 1 , cosh 2 , 



1 



1 _ e -0\k o \ 



(A-8) 



and k is the renormalization point. The "(1,0)1" in V^ 1 ' ) 1 indicates a function 
which has a one-point type 1 external field and a zero-point type 2 external field at 
the one-loop level. —iV^^ e , —iVc and — iVp~ are contributed from Figs.^ 
(a),(b)|| (c)-(e) and Fig|| respectively. 

The above quantities are then integrated over 93, yielding the following: 



V 



1 M 4 ,3 , M 2 . 

f- - In— ) + 



, , , , . rO0 A^2 ln(1 _ -/3 E) (A . 9) 

(4vr) 2 4 v 2 K Jo (2vr) 2 /3 V ; ' V 1 



Vr 



2bubble 



2 * 



AM 2 



2, ,2 



+ 



M 2 M 2 M 2 

^p{a + (7 - In 4vr - 3 + In -^) In 

r(hi- 



M 2 



7T 



2 + -=)J t + J, 



2(2vr) 2V ^3 

2,„2 



AV 



(A10) 
(ATI) 



(AT2) 



where, 



E = \/k 2 + M 2 



M 2 M 2 

(1 - In—) - If 



K 



(4vr) 



(ATS) 
(AT4) 

(ATS) 



2(4vr) 2 

M 2 , , , 1 M 2 x , M 2 , 

_{ d _( 7 _ ln47r _l + ln _ } 

--(7 2 - 2 7 + 2 + — ) + (7 - 1 - - ln4vr) ln4vr ~ -5.68496, (A- 16) 
2 6 2 

7T 2 9 1 

12 + 2 + 2^ 7 ~ ln4 ^ 7 " ln47r " 3 ) 







+ / dx dy — lny(lna 



a 



10.16186 - 2.17195 = 7.98891, 



-(3,-1 + -J— ), /? = -(! + ' 



j; 2 — x' 



a 
It 

U lo (2vr) 4 2T,T 2 
and 7 is the Euler constant. 



x 2 — x' 



dk k 2 n B (E) 
(2T) 2 E ' 

dkidk 2 k 1 k2n B (Ei)n B (E 2 ) 



In 



\E\E\ - (M 2 + 2hk 2 ) 2 



\E\E\ - (M 2 - 2k l k 2 ) 2 



(AT 7) 
(ATS) 

(AT9) 
,(A20) 
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Fig. 




X 

(a) (b) 

One-loop effective potential diagrams. 




(a) 



(b) 



Fig. 7. (a) One-loop diagram contributing to the 1-point function, (b) Counterterm contribution 
(which includes all 7i-order contributions in principle). The index 1 indicates a type 1 vertex. 
The cross represents the counterterms. 








(b) 

Fig. 8. (a) and (b) 2-bubble diagrams, (c)-(e) Tadpole diagrams including counterterms. 




1 1 





(a) 



(b) 



(c) 





Fig. 9. Diagrams which contribute to the setting sun diagram. 



Appendix B 

OPT Equations 



Here, we summarize the calculations of the effective potential and its derivatives 
in the full OPT. 
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(a) (b) (c) (d) 



Fig. 10. Two-loop effective potential diagrams. 




(a) (b) (c) (d) 



Fig. 11. (a) Additional 0(5) contributions to the 1-point function in the full OPT. (b)-(d) Addi- 
tional 0{5 2 ) contributions to the 1-point function. 




(a) (b) 



Fig. 12. (a) Additional O(S) contributions to the effective potential in the full OPT. (b) Additional 
0{5 2 ) contributions. 



The two-loop effective potential in the full optimized perturbation theory reads 
V = V° + V S + V 5 \ (B-l) 
V° = \m\ 2 + | <f\ (B-2) 

irS 1 2 V 4 1 M\3 , M\ 

V =-2™ - & -(4^ — ( 2" ln ^ 

+ r^^-^ (b - 3) 



where, 



2 2 / 2 2\ 2 

\i = m — (m — \x ) = m — Xi 
A = 9 - {g - A) = g - rj, 



M 2 M 2 
^=2W (1 - ln - r) - /i ' (B ' 6) 



(B-5) 
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T f ^S s + C f 

M 2 . n 1 M 2 M 2 

- {c+ (-2 + -ln— )ln— } 



(4vr) 4 ^ ~ 2 ' 

1 , M 2 7T x 



2(2vr) 5 



it 
/i 



f°° dh 
7o (2tt 



d& k 2 n B (E) 



~ Jo ~W 



(2vr) 2 £ ' 
dk\dk 2 kik 2 n B (E 1 )nB(E 2 ) 



1E\E 2 



In 



4£ 2 £| - (M 2 + 2fcifc 2 ) 



4£ 2 £ 2 - (M 2 - 2feiA:2) 2 



(B-7) 

(B-8) 
(B-9) 



M 2 =m 2 + ^ 2 , £?= Vfc 2 + M 2 , n B (£) = (e^-1) -1 , (B-10) 



1 /-i /-i 0--y)P 

c = -(7 — 7 + ln47r) — / cte / lny(lnct ) 

2 Jo io a 

~ 2.30495, 



a = -(y-l + -J—), = -(l + 

x z — X 

and 7 = 0.5772 

B.l. For PMS conditions 

Here, we give equations needed to study the PMS conditions. In the following, 
we set g = A, so rj = in Eq.(B-l): 



x 2 — x 



)> 



(B-ll) 
(B-12) 



dV dV° dV 5 dV 52 



dip dip dip dip 
dV° 
dip 
dV 5 



2,^2 

6 



= ¥ >[- x -A# t ], 



dV 



5 2 



dtp 



w ,, ^ a 2 /rr , a<^ 2 ar f 



dM 2 



2 <9M 2 ' 



The partial derivative with respect to m 2 is found to be 



dV dV° dV s dV 



s 2 



+ 



dm 2 dm 2 dm 2 dm 2 
dV° _ 1 2 
dm 2 ~ I 9 ' 



dm 2 
dV &2 
dm 2 



= K t + ( X + \K t )^ T ± + - 1 - ---- 



5M 2 



4 dM 2 ' 



(B-13) 
(B-14) 
(B-15) 
(B-16) 

(B-17) 
(B-18) 
(B-19) 
(B-20) 
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Also, when <p = 0, (B-17) read 



dV 



dm 2 



r _r [x+XK,) m-- (B ' 21) 



B.2. For FAC conditions 



The second derivative of the effective potential with respect to ip for the FAC 
conditions is found to be 

d*v_cPv°_ an^ d 2 v s ' 2 

dtp 2 dp 2 dp 2 dp 2 ' 
^-m> + ^, (B-23) 



d^ 2 1 * dM 2 ' 

4^ + ^ 2 ^nAV^). (B.25, 



B.3. .ffig/i T forms 

Here, we give the high T {j3M — > 0) forms of Eq.(B-6), derivatives of Eq.(B-6) 
and Eq.(B-9). We find 

T 2 MT M 2 , T 2 M 2 . 1. „„, 

*' - " 24 + 1ST + 2W (1 " " WT'^* ~ 7 + 2 WB - 26) 

8K f T 1 T 2 1 1 

' ' +^2( 1 - ln ^)-77TTo(ln47r- 7 +-), (B-27) 



dM 2 16vrM 2(4vr) 2V k 2 y (4vr) 2 v ' 2 

d 2 ^ T 



8{M 2 ) 2 327rM 3 ' 



(B-28) 



12^)2 0* + 3-48871). (B-29) 
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